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A rigorous derivation of multicomponent diffusion laws 

Zaibao Yang* Wen-An Yongj Yi Zhn^ 


Abstract 

This article is concerned with the dynamics of a mixture of gases. Under the assumption that 
all the gases are isothermal and inviscid, we show that the governing equations have an elegant 
conservation-dissipation structure. With the help of this structure, a multicomponent diffusion law 
is derived mathematically rigorously. This clarifies a long-standing non-uniqueness issue in the field 
for the first time. The multicomponent diffusion law derived here takes the spatial gradient of 
an entropic variable as the thermodynamic forces and satisfies a nonlinear version of the Onsager 
reciprocal relations. 

Keywords. Multicomponent diffusion laws, conservation-dissipation structure, Maxwell itera¬ 
tion, Onsager reciprocal relations 


1 Introduction 

Multicomponent diffusion plays a dominant role in many industrial and natural processes and has 
been extensively studied in a large number of chemical engineering books and articles since 1948. See 
[iiiaiaiiiisiiiTiiiiiiiioiiiiiiiiiiiaiiiiii] and the references cited therein. A fundamental task of 
these studies is to determine the relations between the thermodynamic fluxes and forces, which is called 
multicomponent diffusion laws or constitutive relations. In the literature, various different constitutive 
relations can be found and they are not always consistent with each other. For example, in [9] pub¬ 
lished in 2006 one can find the statement “The fact that Dij (diffusion coefficients) is fundamentally 
nonunique is clearly documented in the literature.. Therefore, there has been a urgent need for 
a mathematical clarification of the situation. 

In this paper, we present a mathematically rigorous derivation of a multicomponent diffusion law 
for a mixture of gases and thereby clarify the uniqueness for the first time. For the sake of simplicity, we 
assume that all the gases are isothermal and inviscid. Our derivation follows the Furry-Williams approach 
[niaiaiio] starting from the macroscopic balance equations of all the gases and uses the Chapman-Enskog 
expansion or Maxwell iteration recently justified in m for a class of hyperbolic relaxation systems. 
The justification is crucially based on a conservation-dissipation structural property of the system of the 
macroscopic equations. The conservation-dissipation structure was firstly proposed and was shown in m 
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to be satisfied by many classical models from mathematical physics. Recently, it has been successfully 
used in [TH] to develop a general theory of mathematical modeling of irreversible processes. 

The multicomponent diffusion law derived here takes the spatial gradient of an entropic variable as 
the thermodynamic forces and satisfies a nonlinear version of the Onsager reciprocal relations. Namely, 
the corresponding diffusion coefficient matrix is symmetric and positive-definite. The use of the entropic 
variable as thermodynamic forces is a key of the original Onsager relation m and provides us a nice 
formulation of the multicomponent diffusion laws. This formulation is mathematically and numerically 
important in treating the multicomponent diffusion systems. 

As mentioned earlier, there are various multicomponent diffusion laws in the literature and it is 
almost impossible to comment on all of them. For most of the existing multicomponent diffusion laws, 
the diffusion coefficient matrix is generally non-symmetric |ini 111) and sometimes even not completely 
determined [5]. A further comparison of ours with that in is given in Section 4. On the other hand, 
the Maxwell iteration was also used in mis] to derive the multicomponent diffusion laws. We also 
notice [Til ITS] where the multicomponent diffusion laws were derived in two different frameworks of 
non-equilibrium thermodynamics. Unfortunately, all these derivations are formal, lack a mathematically 
rigorous justification, and therefore the uniqueness issue has been not resolved before. 

The paper is organized as follows. In Section 2, we introduce the classical mass and momentum equa¬ 
tions for multicomponent diffusion systems. Section 3 focuses on the conservation-dissipation structure 
of the systems. The multicomponent diffusion law is derived in Section 4. 

2 Governing equations 

Consider a mixture of N gases. Assume that all the gases are inviscid and isothermal. Then the 
motion of the gases obeys the classical mass and momentum equations 13 HQ). Let Pi and denote the 
density and mass averaged velocity of species i, respectively. The species-specific mass equations read as 

^ + V ■ (p,v,) = w., ^ = 1,•••,1V, (2.1) 

where Wi represents the net mass production rate (per unit volume) of species i due to chemical reactions. 
Since mass can either be created or destroyed by chemical reactions, the Wi’s always satisfy 

N 

On the other hand, the momentum equation for species i is 

^^^^ + V-(p,V,V,+p,/d)=p.f. + mV, + Gr“, i = lr--,N. (2.2) 

Here pi = Pi{pi) is the pressure of species i, Id is the unit matrix of order d {d = 1, 2, 3), piti is the body 
force acting on species i, and G™** stands for the net impacts of all interspecies collisional momentum 
exchanges on species i. Since interspecies collisions conserve the total momentum, the sum of over 
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all species must be identically zero: 


N 

Qcoiz ^ (-2.3) 

By its definition, is a function ofX=(Vi—Vi,-- - — V^) and vanishes whenever X = 0. 

The latter means that there is no collisional momentum exchange with species i when other species all 
have the same velocity. Based on this observation, we may write 


^coll _ 


ac 


coll 


ax 


-(0X)d0-X = ^Ky(V,- v,) 


with Kij = [K!^f]dxd being d x d-matrixes. is called a collisional coefficient, possibly depending on 
any quantities like densities under consideration. 

For the Stefan-Maxwell model [I1S[9], each Kij is a scalar matrix, that is, 


K 


q:/3 




did — THii Ui' 


Here 5a0 is the Kronecker delta, fhij = rriimj/{rrii +rnj) is the “reduced mass”, rrii is the molecular mass 
of species i, and Vij is the averaged frequency (per unit volume) of collisions between molecules i and j, 
and satisfies 

^ij — O'ji > 0 

for each i ^ j and i, j = 1, ■ ■ ■ , N. 

For future references, we rewrite G™** as 


N 




with 


N 


^ij — ^ij ^ ^ ^ik ' 




for each z, j = 1, • • ■ ^ N. Obviously, Kij = piVi, • • • , Pn'^n) satisfies 


N 

'^Kii = 0 and Kij = Kji <0, i ^ j, (2.4) 

for z, j = 1, • • • , A^. This is consistent with (12.31) . 

In a mixture of gases, collisions usually happen much faster than the macroscopic fluid motions. 
Denote by e the ratio of characteristic collision time to characteristic fluid mechanics time. Then e is 
small and the collision coefhcient Kij may be scaled as 

K,, = (2.5) 


with Kij = 0(1). For the notational convenience, we will use Kij for Kij in the rest of this paper. 
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3 Conservation-dissipation structure 


In this section, we show that the system of equations (EH) and (lO) with fast collisions (1^ possesses 
the conservation-dissipation structure proposed in m To do this, we rewrite (EH) and E21) as 


or 


j=i j 




(3,1) 


1 = 1 


with the external forces ptii and chemical reaction terms W being ignored to simplify the exposition. In 
(EH, d is the space dimension, 


v = 


p^ 

PiVi 


F,(V) = 


Pi^ip 

Pi^iVip +Piej 


( 


Qz(U) = 


0 


N 


- E 


k=l 


u = (Ui,U2,..- F,(U) = (F,(Ui),--- ,f,(u^)E, g(u) = (gi(u),..- ,Qw(u)E, u, 

denotes the j-th component of velocity V^, and stands for the j-th column of the unit matrix I^. 
Following |20] . we define a mathematical entropy rji for species i: 

Vi{pi, Pi^i) = Pi [ 


'Pi 




with Pi a possible positive value. Consequently, a mathematical entropy ? 7 (U) for the mixture can be 
introduced as 

JL JL ^ 

(3.2) 


i=i “'p- ^ i=i 


The conservation-dissipation structure for system EH) can be stated as follows. 

Theorem 3.1 Assume that (12.411 holds true and pi = Pi{pi) is strictly increasing for pt > 0 and for each 
i = 1, • • • ,7V. Then the mathematical entropy iy(U) defined in (13.2|) for system (13.111 is strictly convex in 
Ou = {U S I Pi > 0, z = 1, • • • , N}. Moreover, the following three statements are true: 

d^p{V) 9Fj(U) , ^ ^ . TT ^ 

• —- Trr :— symmetric for each U G Uu and for each j; 

aU 

• There exists a symmetric and semi-positive matrix T(U) such that 


Q(U) = -L(U) 


drm, 

(9U ’ 


• The null-space of L(\J) is independent o/U G Ojj- 

Remark 3.1 The first statement is the well-known entropy condition for hyperbolic conservation laws 
and corresponds to the classical principles of thermodynamics. The second one can be understand as a 
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nonlinearization of the celebrated Onsager reciprocal relation in modern thermodynamics and implies 
the second law of thermodynamics. It displays a direct relation of irreversible processes to the entropy 
change. The last one expresses the fact that physical laws of conservation hold true, no matter what state 
the underlying thermodynamical system is in (eguilibrium, non-equilibrium, and so on). 


Here is our proof of the above theorem. Compute 


9U,; 


/ dVi \ 
dpi 

V V. ) 


d^p{V) 1 
’ dV.dVj ~ J, 


+ -Vf 

-V, h 


Sij. 


d^nCU) 

Because p^{pi) > 0 for pi > 0, one can directly verify that the Hessian is positive definite for 


aU2 


each U S Ojj. Therefore, rj = ry(U) is strictly convex in Ojj. 

Next we compute from the definition of F_, (Ui) that, for i ^ k, 


d'^r](\J) 


dU.dUk 


= 0 , 


gF,-(U.) 

dVt 


= 0 , 


and 


a2,y(U)9F,(U,) _ 1 

dW dV, 


Pi 




d^pfU) dF^(U) . 


Therefore, — " is symmetric. 


9U2 5U 

Now we recall the structure of (5(U) and define a (d+ l)iV x (d+ 1)A^ matrix T(U) = [Tife(U)]ArxAf 


with 


Then (5(U) can be written as 


L^k{U) = 


0 0 
0 Kikid 


Q(U) = -L(U) 


drijV) 

9u ■ 


From (12.4p it is not difficult to see that T(U) is symmetric and nonnegative-definite matrix. 

It remains to show that the null-space of L{\J) is independent of U. To do this, we write w e 

/ Wi \ 


as w = 


with Wfc = 


Ji 


and w" e R<^. If L(U) • 'w = 0, then 


\ Wiv / 


N 


N 


'^Lik'Wk = ^ 


k=l 


k=l 


0 0 
0 Kikid 


wi 


( 


rll 


N 


\ E Kik^i^ 

\ k=l 


= 0 


for each i. Since the null-space of the N x 7V-matrix K = [iFy] is spanned by {(I, • • • , I)^}, it follows 
that 


= • • • = w 
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Therefore, the null-space of T(U) is 


{w e : (x 


\flf) G and 


,11 


= = • ■ 




e 


which is independent of U. This completes the proof. 


4 Multicomponent diffusion laws 

In this section, we derive a multicomponent diffusion law by using the Maxwell iteration [nna 
or Chapman-Enskog expansion justified in our recent paper |16j . The conservation-dissipation structure 
provided a highly efficient framework for the justification. 

To begin with, we introduce the mass diffusion flux (also called thermodynamic flux) 

J, =p,(V,-V), 

where the mass-averaged velocity V is defined as 

N N 

^ i=l i=l 

From these definitions, it follows immediately that 

N 

i=l 

This is called zero-net-flux condition in some literature |2]. Due to this condition, we only need to consider 
the first {N — 1) fluxes 3i (i = 1, • • • , — 1). 

Referring to m, we rewrite the multicomponent diffusion system in term of the new variable 

W = {p,pY,pi, ■ ■ ■ ,pjv_i, Ji, • • • , ^— U = (pi,piVi, • • • , pn, PnYn)^. 


Firstly, the N mass equations are equivalent to 


^+V.(pV) = 0, 

^ + V-(p,V-f J,)=0, t = 


(4.1) 


which is derived by summing up the mass equations (12.111 over i and using the definition 3i = pi (V^ — V). 
Recall that the external forces and reaction sources have been neglected. In order to obtain the equations 
for the total momentum pV and the first {N — 1) fluxes Ji, we recast the left-hand side of the momentum 
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equation (12.211 as 


+ V • (p,V, ® V,) + Vp, 

= ^ + V • (p,v 0 V) + vp, + V • (J, ® V + V ® J, + £i^). 

Summing up these momentum equations over j, we use (1^ and the zero-net-flux condition to obtain 

^^+V-(pV0V)+V^p,-b^V-(:^^-^)=O. (4.2) 


N 


Moreover, we use ^ Kij = 0 in (|2.4I) and rewrite the collision term 
i=i 


N 


Gr“ = --E^b(v,-v) 


N-l 

(i: 

i=i 

N-l 


J, 

K ^ 


Jat \ 


KiN 
Pj Pn 


e ^ Pi PJV 


iV-l 


1 T 1 

= --E^b(- + —E JO 

e ^ Pj PN 


i=i 

Af-l 


= -- E K,,{-Sji + —)3i. 


Thus, the equation for reads as 


dJ, , 9 (p,V) 


dt 


+ 


dt 

N-l 


-(- V • {pi^ 0 "V) -f Vp^ -f V • (V 0 -f 0 V -f 


= -- E K,,{-Sji + —)3i. 


(4.3) 


Consequently, system dSU has been rewritten as (Hlll-dg. 


N 


On the other hand, we recall (12.4|) that Kij = Kji < 0 with i ^ j and ^ Kij = 0. Then 

i=i 

[.f^zj](Ar_i)x(w-i) is a strictly diagonally dominant and symmetric matrix. Therefore, it has an inverse, 
say /C, and the inverse is symmetric and positive definite. Moreover, we set m 


^ ii — ^ii T and Cij — Pj ^ij 

Pj Pn P 
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for z, j = 1, • • • , — 1. Note that 


N-l N-1 

= V Cu<^ij = V {piSu - + —) 

^ ^ P Pj PN 


Z =1 


Z=1 


Pi 


_ . _l_ 


= <5z 


Pn P PPn 
P i Pi 


N-l 

E 

1=1 


Pi 


Pn P PPn 


ip - pn) 


and 


— Sij. 


^ip^V) 

dt 


+ V ■ (p,v (g) V) =p,(Vt + V • VV) + +{p,t + V ■ (P.V))V 


=p,(Vt+V.VV)-(V-JOV 


due to the second line in 6ID. The flux equation 631) can be rewritten as 

N-l T (Vi T 

J. = -e ^ C^klCki [piNt + V • VV) + Vpi - (V • Jz)V + (J,)* + V • (V 0 J, + J; 0 V + 

k,l=l 

AT-l N 

= -e ^ Ofc/Cfcz ( - ^V ^ p, + Vpi + 0(e)). 

k,l=l ^ j=l 

(4.4) 


The second step is due to 621 ) and the Maxwell iteration, using that = 0{e) indicated by the first 
step. 

Furthermore, we deduce from f|4.4D that 


N-l 


N-l 


Jr = 


—e 


= —e 


E c„kJ- Ya _ Zw + + 0(e) 

E PJ)^) + 0(e) 

. . ^ - p Pj p . . Pn ^ 


fc,Z=l 

N-l 


N-l 




£ i: ft.K„( ^(„i„ - a£i)Za _ + o(c 

-- P Pn P P JsJ / 


N-l 


PjPl VpN^ 


fc,Z=l 

N-l 






-£ V c„)c„(pa-^)(^-—) + o(e) 

P Pj PN 


k,l,j=l 

N-l 


-e ^ 0,fc)Cfe,0,(^-^)+0(e2) 

kjjii 


for f = 1, • • ■ , iV — 1. In the last step we have approximated K. = /C(p, pV, pi, • • • , p^v-i, Ji, • • ■ , JAf-i) 
with 1C = /C(p, pV, pi, • • • , pat-i, 0, • • • ,0), that is, the latter is evaluated at the equilibrium. Truncating 














the expansion above, we arrive at the multicomponent diffusion law 


N-l 


3i — £ ^ ^ Dij(u){^ 


i=i 


Vpj _ VpAf 
Pi Pn 


) 


(4.5) 


N-l 

with Dij(u) = ^ CiklCkiCij the multicomponent diffusion coefficients, which depend only on the con- 

k,l=l 

served variable u = {p,pV,pi, ■ ■ ■ ^pm-iY'■ 

As to the relation (145]) . we have the following remarks. 

Remark 4.1 • The diffusion matrix 


D{u) = [Ai(M)](Ar-l)x(Af-l) = [Cij]{N-l)x(N-l)K^[Cij](N-l)x{N-l) 


is symmetric and positive definite, since [Cij]( 7 v_i)x(Ar-i) is invertible, 1C is positive definite and 
they are both symmetric. 




According to the general theory (see, e.g., m), 
spatial gradient of the equilibrium-entropic force. 


the term (—— - Q^n be expressed as the 

Pi PN 

To see this, we recall (ESI) that 


N . 

^ -ti 


2 = 1 
N 






Thus, we have 

=7?(U)|j,=o = + 

For fixed p and V, we compute 


dpi 




VPAT 

Pn 


i = 1,2,--- ,Ar- 1. 


Consequently, the multicomponent diffusion law (|i31) can be rewritten as 


N-l 


Jj = -e ^ D,,j{u)y 




dpi 


Namely, the thermodynamic fluxes Ji are expressed in term of the equilibrium-entropic (thermody- 

namic) forces V— -, with the coefficient matrix D(u) symmetric. This looks like the Onsager 

dpi 

reciprocal relation m, while the definition of the thermodynamic forces is consistent to those in 
the literature \12[ \JW. 


With the diffusion law (03), the system can be approximated formally by the following 
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second-order partial differential equations 


^+V.(pV)=0, 

+ V • (pV 0 V) + V ^ p, = 0, 
i=i 

^ + V • (p,V) = eV • ^ 


Pn 


for 1. Thanks to the remark above, the last equations can be written as 


/ 

Mt -I- V • G{u) = e 

V 



(4.6) 


N 

with G{u) = (pV, pV (8) V -I- ^ PiV, • • • , pAr_iV)^ and ^ = (pi, • • • , pAr_i)^. This is the second- 

i=i 

order partial differential equations (2.5) in [16] with 


= (5jfediag(Od+i,T>(M))pO« (m), 


where Od+i is the zero-matrix of order [d + 1). It is known from |16j that rj^‘^{u) is strictly convex and 
therefore its Hessian ri^^{u) is positive dehnite. Moreover, it was showed in m that the system (021) 
has a nice entropy structure, which is extremely important mathematically and numerically. 

With the above expression of we can simply show the following proposition. 

Proposition 4.1 The system 031) satisfies the isotropy condition (*) of Theorem 2.1 in m- 

Proof. By Lemma 3.2 in m, it suffices to show that the null-space of the symbol matrix 

d d 

B{u,fi)= B^^{u){r]ll{u))~^^jfk= Y Sjkdia.g{0d+i,D{u))^j^k=dia.g{0d+i,D{u))\f\^ 

j,k—l j,k—l 

is independent of u and ^ € i?'^\{0}. This is clear thanks to the positive definiteness of D{u). 


Thanks to this proposition and the conservation-dissipative structure, we use Theorem 2.1 in |16| 
and get the following conclusion. 

Theorem 4.1 Under the conditions of Theorem 3.1, let s > d/2 -|-1 be an integer. Assume that U{e) as 
initial data for PDFs 031) and u{e) for (021) are in for e > 0, satisfy 

||zi(-,e) -u(-,£)|U = O(e^) 

with u{x,e) the conserved mode ofU{x,e), and all the components corresponding to densities have pos¬ 
itive lower bounds. Then there exist e-independent positive constants T* > 0 and K{Tfi) such that the 
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solutions to PDEs dllD and (HU) with the above initial data, denoted by U^{x,t) and Ue{x,t), are in 
C{[0,n],H%R‘‘)) and 

sup ||w®(-,i) - UE(-,i)||s < K{T„)e^, 
te[o,T.] 

where u^{x,t) = pf, • • • is the conserved mode ofU^{x,t). 


The notation used in Theorem 4.1 is standard: For a nonnegative integer s, is the space 

of functions whose distribution derivatives of order < s are all in and we use ||17||s to denote the 
standard norm oi U G . (^([O, T],X) represents the space of continuous functions on [0, T] with values 
in a Banach space X. 

Finally, we give a detailed comparison with the multicomponent diffusion law derived in 

N 

Remark 4.2 In Lam introduced N numbers oji satisfying ^ Wi 7 ^ 0 and replaced the collision co- 

i=l 

efficients Kij with Kij = Kij + iOiPj based on the zero-net-flux condition. By using the invertibility of 
the rank-one modification [Kij] of the singular matrix he derived the following multicomponent 

diffusion law 

N 

3i — £ ^ ^ piDij dj, i — 1, 2, * • * , N, 

i=i 


where [Dij] = p[Kij 
following relation 


— T)' T) ' O' 

1“^ and dj = V— + (—-^)V(lnp) with p = ^ Pi- It is not difficult to verify the 


P 


P 


N 

E 


^ YfEl 

Pj Pn Pj 


^Pn 

Pn 


between our thermodynamic forces and Lam’s. Obviously, Lam’s multicomponent diffusion law is not 
completely determined in general, for the diffusion matrix [Dij] depends on the arbitrary parameters uji. 
Moreover, it is not clear whether [Dij] is symmetric or positive-definite, while so is ours [Dij (u)]. 
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